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We introduce, and propagate wave-packet solutions of, a single qubit system in which geometric 
gauge forces and phases emerge. We investigate under what conditions non-trivial gauge phenomena 
arise, and demonstrate how symmetry breaking is an essential ingredient for realization of the former. 
We illustrate how a "magnetic"-lens, for neutral atoms, can be constructed and find application in 
the manipulation and interferometry of cold atoms. 
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INTRODUCTION 

Gauge symmetry is at the core of our current under- 
standing of how the fundamental constituents of matter 
interact. With the discovery of the geometric phase [I| 
diverse systems, ranging atomic, molecular, optical, con- 
densed matter and nuclear physics, have been identified 
in which gauge phenomena, in addition to those aris- 
ing from fundamental gauge fields, emerge. More re- 
cently, researchers [2] have, via the application of laser 
fields on cold atoms, engineered Hamiltonians that lead 
to effective "magnetic" -like forces on the atoms. This 
advance has great potential in the control and manip- 
ulation of quantum matter 3 . In particular, its ap- 
plication promises the capability to create ensembles 
of neutral atoms that exhibit exotic quantum Hall-like 
behavior [H [5]. 

In this Letter, we introduce a single qubit model that 
possesses non-trivial gauge behavior and whose labora- 
tory realization may offer novel routes to the quantum 
control, and interferometry, of cold atoms. We present 
results of time-dependent calculations for wave-packet 
propagation to illustrate how, and under what condi- 
tions, geometric gauge forces manifest. We show how 
the proposed system mimics that of a charged particle 
scattered by a ferromagnetic medium. We illustrate how 
an effective "magnetic" lens can be engineered and pro- 
pose possible applications. 



THEORY 



Consider the Hamiltonian 

h 2 



H 



2m 



V 2 + H ad {R) 



(1) 



where H a d{R), the adiabatic Hamiltonian describing a 
qubit, is parameterized by the quantum variable R, and 
which can be expressed in the form 

H ad = U(R)H BO U\R). (2) 

A detailed, time-independent, description of such sys- 
tems has been outlined in Ref. [BJ, but here we exploit 
time-dependent methods to enhance and generalize the 
conclusions of that paper. Laboratory realizations of adi- 
abatic Hamiltonians discussed in this Letter could be 



achieved using the techniques discussed in Refs. [3, 4 . 
Because U is a unitary operator the eigenvalues of H a( i 
are solely determined by H B o that, in this study, we re- 
quire to be non-degenerate. We take the eigenstates of 
H BO as our basis and set H BO = Act 3 where A > 
is a constant and 03 is the diagonal Pauli matrix. We 
choose [BJ for U, 

exp(-i cr 3 — y) exp(-icr 2 £l{x)) exp(i a 3 — y) (3) 

where Oi are the Pauli matrices, 4> is a parameter, R = 
(x, y) are 2D Cartesian coordinates, and 



Q(x) = ^ (l + tanh(^x) 
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We seek solutions of ihdt4>(R) = Hi/j(R) which, ex- 
pressed in this basis, take the form of the coupled equa- 
tions 
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where 



V = A cos(2n(x)) 

Via = exp(- i $ y) A sin(2ft(x)). 



(6) 



They may be solved using the split-operator met hod 0, 
and a wave packet at t = to, can be propagated to t = 
t + St, for a small time increment St. Introducing the 
dimensionless units t = t, (£, rj) = (x/L,y/L), where 
L is an arbitrary length scale, we obtain 



iJj(t + St) = UkeUvUre ^{t) 



(7) 



where Uke is a diagonal matrix operator whose elements 
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cos(A St) - i cos(2fi(£)) sin(A St) 
-i exp(i$ 77) sin(2fi(£)) sin(A St) 



-i exp(-i$ 77) sin(2f2(£)) sin(A St) 
cos(A t) + i cos(2Q(£)) sin(A St) 
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FIG. 1: (a) An incident wave in the open channel whose 
potential energy is given by the solid blue line. For 
£ > the transmitted wave is propagated on the 
potential surface given by the red solid line. The brown 
line represent the total collision energy of the system 
and the green dashed line represents the off-diagonal 
coupling between the two potential surfaces, (b) The 
same system now illustrated in the adiabatic picture 
(gauge). The blue line is the BO energy for the open 
channel and the red line represents the BO energy for 
the closed channel. 



Expressed in these units A is a dimensionless param- 
eter as are <£>, (3. Figure ([T]) provides a schematic de- 
scription of the dynamics generated by Hamiltonian (JlJ . 
In Figure ([2^,) we provide a time series contour plot 
of the probability densities |/| 2 and \g\ 2 . At r = 

we place a Gaussian wave-packet tpg = ( ^ J centered 

V 9o J 

£ = —4.0, 77 = with an initial velocity directed along 
the positive £ axis. In this region (17 — > 0) 



Had = 



A 

-A 



and the wave packet evolves as that of a free particle until 
it reaches the interaction region £ ps 0. The wave-packet 
is illustrated by the blue contours in Figure ([2^). The 
initial kinetic energy of the packet was chosen so pene- 
tration of the potential barrier, illustrated in Figure ([I]), 
is prevented. However, the packet can execute a transi- 
tion into the open channel across the barrier. In other 
words, a transition from the g to / channel occurs in the 
region £ w 0. This is illustrated in Figure (|2^) by the 
red contours that represent the wave-packet, probability, 
contours in the / channel. In addition to distortion and 
spreading of the packet there is a noticeable swerve in its 
velocity as it emerges from the interaction region. 

We define the adiabatic amplitudes, 



/ = / cosf2(£) + exp(— i&rj) sin £!(£)# 
g = g cos f2(£) — exp(i$?7) sin f2(£)/ 



and ip 



obeys the following equation 6 
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where A is a non-Abelian, pure, gauge potential. In the 
region £ — > —00, f2 — > and g — > g. Likewise as £ — > 00, 
n — » 7r/2 and g — > f. This behavior is illustrated in Fig- 
ure ^p) where we present a 3D plot for the evolution of 
\g\ 2 . In the adiabatic picture the open channel amplitude 
g evolves in a constant adiabatic potential —A shown in 
Figure ([T|d) . As long as the collision energy is below 
the threshold for excitation into the upper adiabatic, or 
closed, channel the system evolves on a single adiabatic 
surface. Under such conditions the Born-Oppenheimer 
(BO) approximation to the solutions of Eq. ([9| is appro- 
priate. In this approximation, the projection operator 
Pip = g, is applied on Eq. M to get 
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where A = PAP is an Abelian gauge potential 
with non-vanishing curl and b(x) is an induced scalar 
potential [5]. It leads to an effective magnetic induction 

\B(0\ = ^/3$sech 2 (^)cos(|tanh(/30) (11) 

and mimics that incurred on a charged particle that is 
scattered by a ferromagnetic medium. The magnetic in- 
duction is normal to the plane of the page, and is illus- 
trated by the green shaded area in Figure (J2}d) . In Figure 
(J2^i) we also plot, shown by the dashed line, the trajec- 
tory for the solution of the classical equations of motion, 
subjected to a Lorentz force v x B, where \B\ is given 
by Eq. (111. Comparison of the classical path and that 



traced by the centers of the wave-packets shows good 
agreement. The deflection angle suffered by a charged 
particle that is normally incident on a ferromagnetic slab, 
of finite width, with constant magnetic induction B di- 
rected perpendicular to the plane of the page is[B], 
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where <& = f_ d£B(t;) is a flux density, and k the inci- 
dent wave number. 

In Table[l]we tabulate values of the deflection angles 8, 
obtained by calculating the expectation values < £(t) >, 
< 7](t) > for various values of incident, adiabatic, packet 
wave numbers k and <f>. In that table we show the de- 
pendence of 9 on the choice of the energy gap parameter 
A. At lower collision energies, so that k 2 « 2A, we 
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FIG. 2: (a) Blue contours represent the initial g component of the wave packet probability distribution. At time 
the packet executes a transition into the / channel, shown by the red contours, which is energetically open in the 
region £ > 0. Subsequent to the transition the packet evolves as a free particle but with a pronounced swerve in its 
velocity, (b) 3D plot of the adiabatic gauge probability distribution |g(^)| 2 - In both figures the deflection angle has 
the value tan 9 0.59. In these calculations we took k = 12, $ = 6, A = 200, (3 = 2 and 
t = 0, n = 20, r 2 = 30, r 3 = 40. £, t] range between =f27t. 
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FIG. 3: Plot of probability densities \g\ 2 , |/| 2 when 

2 - 2A. 



both channels are open. The threshold energy k t 
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find that Eq. ( 12 1 provides a good approximation for 9. 
As the energy gap 2A is decreased, for a fixed value of 
k, Eq. (12) is less accurate. However, even at threshold 
k 2 s» 2A there is still fairly good agreement between the 
calculated value and that predicted by solutions of Eq. 
(10 1. When k 2 > 2A the excited adiabatic state is open 



and transitions from the adiabatic channel labeled g into 
/ is energetically allowed. In Figure |3^,) we illustrate the 
evolution of the amplitudes |g(i)| 2 and |/(t)| 2 for the col- 
lision energy where 2A/fc 2 = 1/3. The incident packet, 
in the g channel, bifurcates when it reaches the interac- 
tion region. Because there is sufficient kinetic energy, the 
remainder of the initial packet proceeds along the path 
rj = in the region £ > 0. However, a fraction of that 
packet makes a transition into the / channel, and our cal- 
culations show that the angle of the swerve illustrated in 
that figure is in harmony with that obtained at the lower 
collision energies tabulated in Table [I] Therefore there is 
a state-dependent spatial segregation of the initial beam, 
a hallmark of quantum control. In panel (b) of that fig- 
ure we plot these probabilities for energies 2 A/A: 2 = 1/6 



TABLE I: Calculated values for the deflection angle 9 
are tabulated in the third column. The fourth column 
gives the values obtained using Eq. 
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and now find a small, barely noticeable, remnant of the 
packet in the / channel. In the limit k —> oo (or A — » 0) 
Eq. ([5]) allow analytic solutions and g(t) simply evolves 
as that of a free particle. According to definition Eq. ([8] 
) the initial, adiabatic gauge, packet g(t) makes a transi- 
tion into the / channel in the region £ > 0. This "tran- 
sition" is induced by the off-diagonal gauge couplings in 
Eq. ([9]). The "transition" is simply an artifact of the adi- 
abatic gauge (i.e. different definitions for the scattering 
basis in the two asymptotic regions, £ < 0, £ > 0) and 
does constitute a "real" physical transition. In order to 
better understand the behavior illustrated above we re- 
express the unitary operator U that defines the adiabatic 
Hamiltonian. In Ref. [6] we argued that U can always 
be written in the form 



U{R) = V exp(- 
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where V is a path ordering operator, and A is a non- 
Abelian gauge potential. U must be well-defined (i.e. not 
multi- valued) for all R and therefore Eq. ( 13 ) must be 
independent of the path C, or V exp(— i § dR' ■ A) = 1. 
Gauge potentials that satisfy this condition are some- 
times called a pure gauge and typically have vanishing 



4 



curvature everywhere. Because of relation ( 13 1 we con 



elude that A is encoded in the definition of H ac i and 
since [Hbo , A] ^ gauge symmetry is explicitly broken 
by Hbo- Though A is trivial, in the sense of it being 
a pure gauge, quantum evolution selects and is sensitive 
to the projected A = PAP non-trivial connection. In 
the adiabatic picture the gauge potentials are explicit, 
being minimally coupled to the amplitudes. As k — » oo, 
or A — > 0, their presence simply contributes to a mul- 
tichannel, or non-Abelian, phase in the adiabatic ampli- 
tudes that has no physical import. In contrast, at lower 
energies the system behaves as if it has acquired a non- 
integrable phase factor. The effect is most pronounced 
when the excited adiabatic state is closed. 



significantly distorted. At t\ = 30, where the packet de- 
scribes free particle evolution, it assumes the shape of a 
"shark-fin" as shown in that figure. The width, along 
the rj direction, is significantly reduced from its value 
at To. A dramatic consequence of the proposed "mag- 
netic" lensing effect. Such a lens, if realized, could find 
application as an "optical" component in an atom laser. 
In addition, consider two localized but coherent packets 
spatially separated at t = 0. After passing through the 
lens they meet and interfere. Because of different geo- 
metric phase histories the interference pattern depends 
on the "magnetic" flux enclosed by the paths. One can 
therefore anticipate its application as a novel expression 
of atom interferometry. 



APPLICATIONS 



In the discussion up to this point, we have used time- 
dependent methods to validate and extend the conclu- 
sions given in Ref. jS] in which time-independent meth- 
ods allow exact scattering solutions for Hamiltonian ([lj. 
However, time dependent methods can be exploited for 
more complex scattering scenarios. Following an analysis 
similar to that in which Eq. (12 I was derived, we now 



posit the following form for the parameter 
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where 7 is, in general, a complicated function of j3, k, A. 
Here we set it to have the constant value 7 = 1. Using Eq. 



( 14 ) we propagate wave packets for various values of im- 
pact parameter. In Figure Ha) we plot trajectories of the 
total expectation values < i/}(t)\£\ip(t) >,< ip(t)\r]\tp(t) > 
for the various impact parameters b. At each impact pa- 
rameter we choose identical wave-packet widths and set 
k = 12, A = 400. The trajectories shown in that fig- 
ure, by the solid red lines, demonstrate that the paths 
converge to a common focal point given by / = 3. This 
result is gauge invariant, i.e. it can be obtained using am- 
plitudes obtained in either diabatic or adiabatic gauges. 
However, the adiabatic picture provides a transparent 
physical description. For, in it, the system is accurately 
described by Eq. (10 1. That description includes the 
emergence of an effective magnetic induction B = V x A 
whose magnitude is 



7173*77 (8/ 2 +7? 2 ) sech 2 (/?£) cos (!?Ttanh(/3£)) 
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and is normal to the plane of the page. In Figure (|4J: 
we illustrate the propagation of a coherent wave packet 
slab of finite width along the 77 direction. After its pas- 
sage through the "magnetic" lens at £ ps 0, its shape is 
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FIG. 4: (a) Trajectories of wave-packet expectation 
values for various values of impact parameter. The 



shaded region is a density plot of Eq. (151 for the 
"magnetic" induction, (b) A wave packet slab having 
width d — 2 along the rj axis at To — is propagated to 
the position shown at t\ = 30. 
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